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Replica field theory is used to study the n-dependent free energy of the Ising spin glass in 
a first order perturbative treatment. Large sample-to-sample deviations of the free energy 
from its quenched average prove to be Gaussian, independently of the special structure of 
the order parameter. The free energy difference between the replica symmetric and (infinite 
level) replica symmetry broken phases is studied in details: the line n(T) where it is zero 
coincides with the Almeida- Thouless line for d > 8. The dimensional domain 6 < d < 8 is 
more complicated, and several scenarios are possible. 
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I. INTRODUCTION 

Large scale free energy deviations F = N f of a quenched disordered system from its selfaveraged 
mean Nft yp are extremely rare events with probability ~ e N with the system size N — > do, and 
L(f) = O(l) < for / = O(l) ^ /t yp . Only the case / < / typ is considered in this paper. As invented 
in Ref. [1], L(f) can be computed by the replica method; namely it is the Legendre transform of the 
n-dependent free energy $ of ([I]): 

0- 1 L(f) = -n<5> + nf and -^-(n^>) = f. 

an 

The replica number n is considered a nonnegative real number corresponding to / < /typ, see also 

as 

. In most applications, we have $ — /t yp ~ —n a for n — > + with some positive integer exponent 
a characteristic to the model, and also to the region of control parameters (temperature, magnetic field, 
etc.) considered^ Legendre transforming then provides L(f) ~ —(/typ — f)~°~- In the most common 
a = 1 and the large scale free energy fluctuations are Gaussian. 



situations 



A very atypical behaviour was found many years ago by Kondor 4] for the truncated version of the 
Sherrington-Kirkpatrick (SK) model [5|, i.e. for the mean field Ising spin glass, just below the transition 
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temperature. The exponent a = 5 was later proved to be true in the whole spin glass phase, including 
zero temperature j^, Q]. The original motivation of the present work was to find out how the large free 
energy deviations are influenced by including the geometry of a d-dimensional hypercubic lattice. For that 
purpose, we followed the usual steps to generate a field theory convinient for a perturbative calculation 
whose zeroth order is the mean field (SK) theory. These steps are: 

• Transformation of the lattice system to a field theory, see Refs. [(], Q]. The natural ultraviolet 
cutoff is then the boundary of the Brillouin zone, and it is A = O(l) when length is measured in 
the units of the lattice spacing. 

• The momentum space is shrinking such that T c -T~t<A<1. 

• Momentum dependence is truncated up to the quadratic term. 

We are led then to the replicated field theory presented in the next section, Eqs. ([2]) and ([3]). (Zero 
external magnetic field is assumed throughout the present paper.) It turns out from the calculation of 
Sec. IIVI that the leading behaviour of the large scale fluctuations changes to Gaussian when first order 



corrections to the low temperature SK model are taken into account (this has been noticed in Ref. [§]), 
and the anomalous n 5 term is now subleading. It becomes also clear that this Gaussian behaviour follows 
for any structure of the order parameter as long as replica equivalence is assumed. 

The glassy phase of the SK model has the ultrametric structure with infinite step replica symmetry 
breaking [9| (RSB) 1 . For T < T c , this RSB state is reached by a phase transition 4j when n is lowered 
from around 1 — where $(n) is the free energy of the stable replica symmetric (RS) phase with a nonzero 
order parameter — at some nxr{T) having two characteristics: 



id ] instability, i.e. the so called 



1. The RS phase becomes unstable due to the Almeida-Thouless 
replicon mass is zero for n = ri£r(T). 

2. = ^rsb at the transition. 

Whether or not this feature of the SK model is perturbatively stable is studied in the second part of 
the paper, and rather nontrivial calculations lead us to a positive conclusion for d > 8. The transitional 
domain of dimensions above the upper critical dimension (6 < d < 8) is, however, more complicated: it 
is argued that the question cannot be settled by this first order calculation. 



1 The acronym RSB is used throughout the paper for the special type of replica symmetry breaking with infinite number 
of levels in the ultrametric hierarchy. 
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The outline of the paper is as follows: Sec. [II] presents the replicated field theoretical model suitable 
for the perturbative calculation, whose basic formulae are provided in Sec. IIIIl The one-loop results for 
the re-dependent free energy are shown in Sec. IIVI up to fourth order in the double series of r ~ T c — T 
and n: these are the leading contributions when d > 8. Nonanalytic temperature dependences are more 
and more important while approaching dimension 8, and even for 6 < d < 8. They are computed for 
both the RSB and RS schemes in Sec. |V1 The line of equal free energies of the two schemes is calculated 
for d > 8 in Sec. IVIl whereas we analyze the situation in the transitional regime just above the upper 
critical dimension, i.e. 6 < d < 8, in Sec. IVIli One-loop results for the RS replicon mass Tr are presented 
and the AT instability line deduced in Sec. IVIIIl Some concluding remarks are left to Sec. IIXI Three 
appendices contain some computational details. 



II. FIELD THEORETIC MODEL FOR THE CALCULATION OF THE FREE ENERGY 



In this section, we will define the replicated field theory which is appropriate for the calculation of 
perturbative corrections to the re-dependent free energy $(n,/3) defined as 

1 



nN 



In Z n 
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where Z is the partition function of an Ising spin glass defined on the <i-dimensional hypercubic lattice 
consisting of N Ising spins, and the bar denotes averaging over the quenched Gaussian disorder with 
zero mean and variance J 2 . The thermodynamic limit N — > oo is always understood, whereas re > 
is small but kept finite. We follow the strategy suggested in Ref. ll|: by conviniently choosing the 
bare parameters, the field theory provides the SK results in the tree approximation, i.e. when 

neglecting loops. In zero external magnetic field we can express Z n in the following functional integral 



form (see Refs. 



Z n = C x / [d<t>] e 



(2) 



with the normalization factor C to be fixed later by choosing it to match the SK result. The Lagrangian 
C has the form 



" V ' V P1P2P3 a/37 



1 



Ui 



p ' ap 



24N 



Ypi <Tp 2 ^P3^P4 

a/375 



P2 ^P3 ^P4 



P1P2P3P4 



+«*£ «*?J-- 



a/3 



a/37 



a/37<5 



(3) 
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where the ellipses are for higher order replica s ym metric invariants which are consistent with the extra 
symmetry of the zero magnetic field subspace (id . fl^ ]. In this n(n — l)/2 component field theory the 
fluctuating fields are symmetric in the replica indices with zero diagonals: = </>p and = 0, 
a,P = 1, ... ,n. (Momentum conservation is indicated by the primed summations.) The bare mass m 
depends on temperature as m = — ^j(T' c mf2 — T 2 ) where k is Boltzmann's constant, and /cT c mf = J - 
the mean field critical temperature — differs from the exact one, T c . Introducing the new parameter 
r = 7[ji{Tc — T 2 ) measuring the distance from criticality, the bare mass can be written as m = m c — r 
with m c being obviously one-loop order. Beside r, the couplings u>, u\, U2, 113, U4, . . .parametrize the 
Lagrangian. Around the critical point they can be considered as constants; a complete matching with 
the SK results is achieved by choosing w = 1, ui = 3, U2 = 2, 113 = —6, ua = [11]. 



The spin glass phase below T c is characterized by the nonzero order parameter (f> al3 = -p=(^p£ ) where 
the average is now taken by the measure proportional to e~'~ . It is useful to redefine the fields by the shift 
4>p^ — > (/>p^ — (f) al3 SpL . By this transformation, the new fields continue fluctuating around zero; 
on the other hand, however, the Lagrangian has lost the higher symmetry of the paramagnetic phase 
resulting in the following generic theory (equally convinient for an RS or RSB ansatz): 2 

£ = £ (0) +£ (l) +£ (2) +£ (3) +£ (4) + ____ 



For the sake of simplifying the notations, we will keep the symbols (ftp 13 and C for the transformed quantities. 



The terms above can be worked out using the results of the Appendices B and D of Ref. [13(] providing 



a/3 a/37 a/fyl 



a/37 



a/37 

2 +«4E^V* 2 

a/376 

07 



5a 



af3~{8 
+ ..., 



+« 2 E^ a/34 

a/3 



(4) 



-«4 



A«/3 



76 



+ ... 



(5) 



,2 1 



^ (2) = ^E Ef^+^-^^-^E 

P L a/3 V 



2 1 



7<5 ' 



+ E f - ^ - 5«i E ^V' - ^3^ 7 7/3 + . . . ) x ^5 



a/37 



+ E ( - - ^ a V + • • • ) >< 

a/375 ^ ' 

-s E' fE( to +-) x <^£^ + E(^+-) x * 



(6) 



P1P2P3 



a/37 



a/3 



P2 ^P3 



a/37 



+ E + • • • x + E +•■■)* <^p^ p . 



a/37<5 



+ EK + - 



k a P AP-l rhP s 



P2 ^P3 



a/37<5 



+ ... 



(7) 



£( 4 ) has been omitted here, as it keeps its form in Eq. ([3|) up to this order, whereas the ellipsis dots 
indicate the terms higher than quartic. 

III. PERTURB ATIVE CALCULATION OF THE FREE ENERGY 



Two steps lead to a systematic perturbative treatment: 

• Firstly an interaction Lagrangian is detached, and £^ identified as the non-interractive part: 

£ = £(0) + ^E E (/fc + ^)^ P +^ 

P (a/3), (7,5) 



(8) 



where in the second term the bare mass is now represented by the mass operator M a g~s, and the 
summation Yl( a j3) = Sa</3 * s over * ne n i n ~ l)/2 pairs of replica indices, i.e., over the independent 
field components. Comparing with Eq. ©, the three different types of mass components are simply 
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derived: 



2 



1 1 

M a p, a p = 2m - + 3u 2 + 2u 3 + 4u 4 ) ^ - -{2u x + u 3 + nu 4 ) V ^ ' + ■ 



M, 



-w <j><# - \ («i + u 3 + 4^ 4 ) a7 / 7 -kVfV { + -, (9) 

M Q/3i75 = (^V 5 + 0<V*) - ^4 ^V 5 + • • • • 

Replica equivalence is a must even when replica symmetry is broken , and it has been assumed 
in the above formulas; see the last term of the diagonal component which is in fact independent 
of (5. The free propagator of this generic RSB field theory is the n(n — l)/2 x n(n — l)/2 matrix 
G = (p 2 + M)- 1 . 

The order parameter (p a/3 must satisfy the "equation of state" (</3p£ ) = 0> U P to one-loop order it 
takes the form 

H «P + \ ^.7^ x Jj ^2 G-y6,nv(p) + 2-loop terms = 0. (10) 

(7<S),(A«0 P 

In perturbative computations, the capital letter representation of the couplings is more natural, 
the one- and three-point ones can be deduced by comparing with Eqs. © and ([7]): 



2 6 V 

7 7<5 7 

c2 



H a p = 2 X 



and the eight different types of cubic couplings: 



W aj3 ^ aa =W + W a p )a p )0l p = 2(ui + U 2 + U 3 + 2u A )(j) a P + 

W a p, a p,p^ = i(2ui + u 3 + 4n 4 )/ 7 + . . . , W a p >a p„ s = + . . . , (12) 

W a p )0rft ps = ^u 1 4>' y5 + ... 

where the dots are for the next — 0{4> 2 ) — order, and the three missing cubic vertices — W a /3 :a ~ /ja s, 
Py,ixu, an d W a p^s,iJ,u enter only at this higher order. 

The critical temperature of the field theory, or equivalently m c , can be simply deduced from Eq. (I10p 
by expanding it up to its leading 0{4>) order after setting r = 0. Making use of the large momentum 
expansion of the free propagator G, together with Eqs. ([9]) and (|12p . it straightforwardly follows: 

mc = \{n-2)w 2 ^- \ + l\(2n-l)u 1 +3u2 + {n+l)u3 + (n 2 -n+4)u 4 \^- T V^ + 2-loop terms. (13) 
2 TV ^-^ r> N v 

p p 
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The free energy lnZ n has the following expansion when C l is handled as a perturbation; see ([2]) and ([8]): 

lnZ^ = In C - + In Z G - {L l ) G + \ ((L 1 2 > G - (L%) + ..., (14) 

and the Gaussian averages {■ ■ - )g are taken by the measure e~ c< ^ /Zq- As "tadpole" diagrams are missing 
now due to Eq. (110|) . the only one-loop term is the Gaussian free energy lnZ G = J[d(j)]e~ c<2) which has 
the familiar form in terms of the eigenvalues Ay's of the mass operator M: 

n(n-l)/2 2 r 

l^c = -±E E (15) 

P 3=1 

Using the following identities: 

E^+^^MEO+^E 1 



p 2 ; V a 2 / ; J d ^ p 2 + A 

p p p 7 



P 2 + Aj p 2 p 4 p 6 (p 2 + Aj) 

where A is the ultraviolet cutoff, and frequently applying the trivial relationship Y^ p 1 = A k Y^ p ^ 
with d > k, the Gaussian free energy can be arranged into its final form for d > 8 

^(H>4)*E'-to)*Ei+te)*Ei 



d 2 vr / ^ 2\^i J J ^ p 2 3 J ^ p 4 

P 3 P 3 P 

54 (16) 

-KE\ 3 )xEi + 5(E^^xEi^EE^TT+o(\ 5 )xEi- 



3 P 3 P 3 P V 

Before displaying the result for the free energy, it must be noticed that the contribution Nn^m c ( 



i2\aa 



from £(°) 3 is exactly cancelled by corresponding terms of InZc] see Eqs. @, (|13p and the second terms 
of the right hand sides of (|Alj) . (|A2j) . Moreover, as is stationary at the zero-loop order parameter, 
we do not need to compute the one- loop correction to 4> ■ Substituting the traces from App.|A]into Eq. 



Matrix notations, like ((f> 2 ) aa here, are frequently used in the following part of the paper. 



s 



16]) . the replicated free energy in (|14p takes the following form: 
1 . 1 



l para 



«/3 4 



\nZ n = \lnZ 

nN nN L 

+ ir x (0 2 ) Qa + _ 1 [ui + 3u 2 + (n - l)u 3 + (n 2 - n - 4)u 4 ] I 4 - 2(n - 2) W 2 / 6 } x (</> 3 ) aQ 

+ 7^{ n i ~~ — 3)ui + 6^2 + 2(n — 1)«3 + 2(re 2 — n — 8)1x4] I4 — 2w 2 [(2n — 7)u% — 3u2 — (re — l)us 

- (n 2 - n + 4)1*4] ^6 + 3(3n - 4)w 4 / 8 } x (</> 4 ) QQ + — |u 2 + - [2u 2 - 4(n - 2)uiu 2 + 3u| - 2(n - 5) 
- 2(n 2 - n - 8)u 2 u 4 + 2u|] I 4 - 4w 2 [ - 2ui + (n + l)u 2 + 2u 3 ] J 6 + 24u> 4 / 8 |x^] 

+ ^{( n 3 + n U4) + - [(3n — 4)u 2 + 12uiU2 + 8uiu^ — (n — 3)ii 2 — 2(n 2 — n — 8)^3^4 — n(n 2 — n — 8)n|]/4 
- 2u> 2 [ - (6n - 19)m - 3u 2 + (n - 5)u 3 + (n 2 - n + 4)u 4 ]l 6 + 3(n - 12)u> 4 / 8 } x ((<j) 2 ) aa ) 2 . (17) 
Some remarks are appropriate here: 

• Only replica equivalence was used to derive the above formula which gives the first order — i.e., 
one-loop — correction to the mean field free energy. Any special RSB scheme is included in the 
"invariants" like ((p 3 ) aa . 

• For = 0, the analytic continuation of the paramagnetic free energy is obtained: 



A 



[ln ^para = lnC7 _n(n_l) ^ ln P - 2r 



7T 

p 

an expansion of this formula gives the r terms in InZc, see f)16j) and App.(XJ n = 1 is the annealed 
model; the disappearance of the loop corrections expresses the triviality of this case: In Z = 
iVln2 + ^B. 2 ^y 2 , where the number of interactions Ng can be expressed by the coordination 
number z as Ng = N%. This gives for InC: 



1 J 2 
lnC = ln2 + ' 



nN *~ ~ ~" ' 4(fcT) 2 
which agrees with the correct value of the SK model for n generic. 

The notation Ij. = Yip ^ was used in (117j) which is a truncated formula in the following sense: 

1. The nonanalytic term proportional to r^/ 2 , see the last but one loop integral in (|16p . was 
neglected, as it is subleading in the region of space dimensions studied here, i.e. d > 8. 

2. Higher than quartic terms in (|17p were left out for the same reason. 
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3. We restrict ourself to the model where quartic couplings are the highest order, i.e. neglecting 
the invariants in ([3]) which are there represented by the dots. A fifth order coupling, for 
instance, would provide a contribution ~ I2 x (0 3 ) oa in Eq. (fl~7l) . 



IV. n-DEPENDENT FREE ENERGY IN SPECIAL CASES 

Eq. CEI]) is now used to compute /3A$ = - ln ^™ ~ [hi^] pam |, i.e. the shift of the n- 
dependent free energy from the continuation of the paramagnetic one, in some special cases. 



A. Neglecting loops: the tree approximation 

Inserting the expansions from (|B5|) into (|17|) . and neglecting loop terms, the mean field free energy of 
the model with nonzero r, w, u±, 112, U3, U4, and assuming infinite step RSB takes the form: 

/3Ad> mf = i«rV + + i«a - ^urV + ^ " ^^" 4 nr 4 - ^ufw±n\ (18) 

Inserting u; = 1, ui = 3, U2 = 2, U3 = — 6, and U4 = 0, this gives the free energy of the SK model 

+ v + ii T 4 +0(T 5 ) _^_„5 > 2T= i_f^y. 

6 24 v 7 5120 V J / 

This formula agrees with the results of Ref. [3] up to the order studied here. 4 From Eq. (|18p follows that 
infinite step RSB is only a necessary condition for the anomalous (non-Gaussian) free energy fluctuation: a 
nonzero 114 produces a term linear in n. This may be a generic feature: interactions which are disconnected 
in replica space (like Y^afHS 4> a ^ 2 <P s2 ) generate Gaussian free energy fluctuations even at the level of the 
tree approximation. 



(3$ 



SK 



In 2 + 



4(kTY< 



B. One-loop correction for the infinite step RSB case 

When substituting U3 by u% = U3 + nu^, the n-dependent free energy shift takes the following simple 
form in the first order perturbative calculation [see Eqs. (fl~T)h (|B5p . (fl"8j) ]: 

/?A$ = /?A$(°) + n /3A$W + /?A$ anom , 



4 Unfortunately the definition of r differs from that of Refs. 0, [3, 4], which is called here r'. The simple relation between 
them is r = r' — \t' 2 . 



with 
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/3A$ 



(°) = -urV \l + (ui + 3u 2 -u 3 - 4u 4 U 4 - 8w 2 I 6 ] + -^w~ 4 t 4 





24 



(3ui + u 2 - u 3 ) 



1 



(33uf + 38uiU2 — 26u\u 3 — 24u\u 4 + 21u 2 — 14u 2 u 3 — 8u 2 U4 + 5u| + 8u3U 4 )I 4 



+ 8(5uj - u 2 - u 3 + 4u 4 ) w 2 / 6 + 24w 4 / 8 , (19) 

1 1 r 1 

/3A$ (1) = -tiTV fc 3 I 4 + 4w 2 / 6 l + — a w~ 4 t a -(-9u 2 - 4niu 2 + !2 Ul u 3 + 8nm 4 + 4u 2 u 3 
24 

- 8u 2 u 4 - 5u§ + 81$ J 4 + 4(ti 2 + 4u 4 ) w 2 / 6 + 24w 4 / 8 ] , (20) 



and 



/?A$ anom = -^^ 2 4 w 4 {^2 + [ - 2n 2 + 2(2n - 5)uiu 2 - 9U 2 , - 8u 2 u 3 - 8u 2 u 4 - 2u 2 3 ]h 



+ 4[ — 6«i + (n + 7)u 2 + 6u 3 — 6nu 4 ] u> 2 /6 — 72 w 4 /g| 



x n 



(21) 



C. For comparison: the replica symmetric free energy 

The RS mean held order parameter q = (j) a ^, for any a 7^ (3, satishes Eq. (|B1|) which now takes the 
form 



t = — ~(n — 2) wq — ~ T(n 2 - 3n + 3)«i + u 2 + (n - l)u 3 l t/ 2 . 
2 b 



This equation can be used to substitute r for g in the RS relations (</> 2 ) QQ: = (n — l)g 2 , (</> 3 ) a ° = 
(n — l)(n — 2) q 3 , (</> 4 )° Q = (n — l)(n 2 — 3n + 3) q 4 , a ^ 4 = (n — 1) g 4 . It is now straightforward to 
derive from Eq. (|17p the free energy shift of the replica symmetric system with respect to the paramagnet: 



/3A$ RS = /3A$gg + n /3A$« + n 2 0A$™ + .... 
We can make the following observations: 

(k) 

• The RS free energy is a regular power series in n, with no anomalous part, and the terms /3A$pg 
all have the same character as a power series of r starting with r 3 . 



0A$ { S = /?A$(°) of Eq. (JT2D up to r 4 (inclusively). 



The leading term proportional to r 3 of /3A<I>^g is identical with the corresponding RSB contribution 
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in Eq. (|2DJ). We have 

= ^w~ 2 t 3 [u 3 I 4 + Aw 2 I b ] + ^w~ 4 t 4 \u 2 + \{~7u\ + W Ul u 2 + !2 Ul u 3 + 8niu 4 
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+ 2lu\ + 8u 2 u 3 - I6u 2 u 4 - 3«g + 8uf) h + 2(4ui - 8u 2 - 4u 3 + 8u 4 ) w 2 /e + 48w 4 / 8 

• From the previous observations follows that the RS and RSB free energies differ only at the 0(5) 
order in the double series of r and n. 

The 0(n 2 ) term — which is missing in the RSB scheme at one- loop level — is given by 

(3A^l = ~^ w ~ 2t3 [1 + (ui + 3 ^2 - u 3 - 4n 4 ) h - 8 w 2 I G ] + O (w~ 4 r 4 ) . (22) 

V. THE NONANALYTIC TEMPERATURE DEPENDENCE OF THE FREE ENERGY 

As we have seen in the previous section, the RSB free energy starts to differ from the RS one only in 
the fifth analytic order 0(r 5 ), 0(r 4 n), 0(r 3 n 2 ), and — lastly but most importantly — in the anomalous 
term 0(n 5 ). In dimensions not much higher than 8, however, these terms are subdominate with respect to 
the nonanalytic one proportional to T d l 2 . For retrieving this nonanalytic contribution from the Gaussian 
free energy in Eq. (|15p . it is converted into the following equivalent form [Eq. (|16p is not very useful for 
that purpose]: 

,2 



P 



3 P J P 

and the last term gives the nonanalytic contribution In which can be generally written as 

— lnZ?, a = - — G aBaB . (23) 

P (a0) 

A. lnZ£ a for the infinite step RSB sheme 

At first sight, computing the right hand side of Eq. (|23|) seems to be a formidable task due to the 
complicated structure of the Gaussian propagators; see Ref. jl5] for details. Nevertheless, as we are 
interested in the leading nonanalytic temperature dependence, some important simplifications follow: 

• This leading term comes from the near infrared region p 2 ~ r; the propagators there were all listed 



in Sec. 6 of Ref. 



When computing this near infrared propagators, the quartic couplings can 



be neglected in the masses in Eq. Q, they enter only through the order parameter 
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The n = propagators — which were computed in [15j — may be used, as their n-dependence 
enters only at much higher orders. [This is just alike the order parameter function q(x).] 

Beside the replica summations in Eq. 1231 the origin of the relevant n-dependence is xq; see (|B4p . 



The trace in (|23h can be written as an integral over the continuous overlap parameter x: 

(a/3) G a /3,al3 



n ^a/3,a/3 = — \ Gff dx and, by using Eq. (60) from Ref. [15], we can write 



/„' " 2 <5S * - a + * + 2f)(i - ») - a i±|±* ^ J_ £ q(xf ix + _J_ w > J_ £ q{xf dx 

where the dimensionless variable t = wqi/p 2 was introduced. Eq. (jB5[) can be used to compute 
f^q(x) 2 dx = —J2/3 c l )a ^ 2 = — (0 2 ) aQ > and f^q(x) 4 dx = —J2/3 c l )a ^ 4 - Keeping only the relevant con- 
tributions, 5 we arrive at 

£ " 2 * - 1 + 2t - 843 + 32 £w - (» + a + 2(2 > " + s gk? *" V:4 " 5 ' 

We are now in the position to detach the leading nonanalytic temperature dependence in Eq. (|23f) : 
approximating wqi ~ r, 

for d > 8 infinite step RSB : 

(24) 



oo 

^l nZ na__16 / P 2 + l x ^/2_^_ 7 dtp 1 « -4-/2-4* 

niV _ d7 (27r) d p6 (;p 2 + 2)2 xr 20d _/ ( 27r )dp4( p 2 + 2 )2 xwn 2 



[Dimensional regularization is to be understood here, allowing the ultraviolet cutoff A go to infinity. This 
same remark is applicable for the RS case, Eq. (|27|).] Note the lack of the r rf / 2 n term in the RSB scheme. 
The 0(r rf//2_4 n 5 ) contribution is negligible; see Eq. ([2"T]) . 



B. In Zg a for the RS case 



Eq. (|23|) is now simplified as 



nN 2d N 

p 



with the diagonal propagator G\ satisfying 



1 , . - l l l . , l 

-n(n - 1) Gi = - + [n - 1) - + -n(n - 3) - 

2 + A L + Aa 2 + A R 



5 Terms with an extra xi factor are subleading, and thus unimportant here. Their calculation using the near infrared 
propagators from Ref. [3] would even be inconsistent, since corrections to these propagators give the same (subleading) 
order. 
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where the three masses are in leading order 

4 - 2 

A L = 2r, = -t, and A/j = -n -r. (26) 

n — 2 n — 2 

Extracting the terms providing nonanalytic temperature dependence, Eq. (|25[) takes the following form: 

for d > 8 RS : 

oo oo 

(27) 

J_l nZ na = _16 7 d*p p 2 + l rf/2 8 ? P 2 + l d/2 

niV G " d 7 (2^) d p 6 (p 2 + 2) 2 dy (27r) d p\p 2 + 2) 3 

A comparison of Eqs. (|24p and (|2T[) shows that the leading nonanalytic free energy contributions of the 
RSB and RS phases coincide, and the free energy difference is 0(5) in the double series of r and n down to 
dimensions 8. In the next section, the line in the r, n plane where the free energy difference disappeares 
is computed in first order perturbation theory taking into account the 0(5) analytical terms. 

VI. THE LINE OF EQUAL FREE ENERGIES OF THE RSB AND RS SCHEMES FOR d > 8 

To find the fifth order (in the double series of r and n) results, we must do two further steps: 

• The one-replica quantities in Eq. (fT?]) must be extended to the appropriate order by using the 
formulae in appendix [Bl (for RSB) and in subsection IIV CI (for RS). 

• Eq. (|17p must be supplemented by the fifth order terms. 

As for the first step, the notation 5 . . . is introduced to mean the difference of a one-replica quantity 
in the two schemes; 5{4> 2 ) aa = (</> 2 )rsb — ( ( / )2 )rs f° r instance. Then we have 

2) «« = _I W -6 U 2 r 4 + 1 -4 U2 nT 3 _ 1-2 ^2 + Q ^ 

9 3 4 

2 1 27 

3 ) m = ^w~ 7 u 2 2 r 5 - w~ 5 u 2 nr 4 + -w~ 3 n 2 r 3 + -^w 3 u 2 3 n 5 + 0(6), (28) 

5 y ^ = -^-w- Q u 2 r 5 + w~ 4 nr 4 - ^w% 4 n 5 + 0(6); 
— ' 15 80 

whereas 8{4> 4 ) aa and ((</> 2 ) aQ ) 2 are of order 0(6). 

Inclusion of the three fifth order one-replica quantities (cj) 5 ) aa , ^pfp"^ {4> 2 ) a ^ \ and (4> 2 ) aa (<f> 3 ) aa into 
Eq. (|17p is not an easy task. Fortunately, however, the <5's formed of them are at most of order 0(6), 
and thus negligible for the present purpose. 
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We can now search for the line in the r-n plane where the RS and infinite step RSB free energies 
coincide. Inserting n = cw^ 2 u 2 t into Eq. (I28p . the following — somewhat surprising — result turns up: 

6(<j> 2 ) aa = -\(c-2/3) 2 w- e u 2 2 r\ 

8(cf> 3 ) aa = JL( C _ 2 /3) 2 (3c 3 + 4c 2 + 4c + 8) w~ 7 u 2 2 r 5 , 

o(J 

5^2 <p a ^ = - A( c _ 2/3) 2 (27c 3 + 36c 2 + 36c + 32) u;- 6 u 2 r 5 . 

/3 

We can write the free energy difference between the RSB and RS phases up to one-loop order, using Eqs. 
([1]) and (|17p . in terms of c as follows: 



c - ~ - 2u 2 - 1 /d(A) 



where the correction term is 



f d (A) = -(2u\ + lluiu 2 + 12u| + 7u 2 n 3 + 4u 2 u 4 + 2u§) I 4 + 4(2«i - 3u 2 - 2u 3 )w 2 / 6 + 24w 4 / 8 . (30) 
3 

It can be seen from Eq. f|29|) that for the cases c ~ 2/3, the RSB and RS free energies differ only at 3- loop 
order! Without any further assumption, we should go up to that order to gain the 1-loop correction to 
c eq , i.e. to the c value of the line where the two free energies are equal. We will, however, assume that 
for d > 8 the mean field form of Eq. (|29p . i.e. 



j 9(* HSB -* Ra )~(c- Caq ) 3 , (31) 

remains valid for the model. Eq. (|29|) then provides c eq = |[1 + u 2 ~ 1 /d(A)], and finally 

for d > 8 $ RSB = $ RS => n = - [u 2 + f d (A) + 0(2-loop)] w~ 2 t. (32) 

3 

VII. BETWEEN SIX AND EIGHT DIMENSIONS 
A. The infinite step RSB case 



As it has been explained in Ref. [151 ] . the quartic vertex whose bare value is w 2 6 suffers a change in 
its temperature scaling from t° to t '/ 2-4 when crossing d = 8. The RSB order parameter is severely 
influenced by this, as best seen from Eqs. (I10p . (jlip by extracting a term proportional to q(x) 3 from the 
one-loop contribution, see and matching it to the corresponding zero-loop one: 



d d p 1 



3L J (2vr) d p^{p 2 + 2) 2 



x 6 a ^. 



In [2] u = U2/2 was used. 
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The second part between the brackets is called U2 and reproduced in Appendix [C] it obviously dominates 
the bare value Ui- Beside u 2 , the combination u\ — |tt3 emerges too from the one-loop correction, as 
shown in Appendix O and the equation of state in (jlOp and (|11|) . with the relevant terms kept only, takes 
the following form [the usual infinite step ultrametric ansatz is assumed here, and therefore q(x) replaces 



2rq(x) + w 



2 

n 

3 



^o)<?o - 2 (<?i - ^x 1 qxjq{x) - ^xq{xf 



+ { ui - ^ 3 j qh(x) + \u 2 q(xf = 0. (33) 



Solving this equation provides the order parameter function in this transitional domain of dimensions: 

q(x) = ^-x, T = wqx-- (in + u 2 - -U3) q 2 , n = -xq, 6 < d < 8; (34) 
U2 2, \ 3 J 3 

compare it with (|B2j) . (|B3p . and (|B4|) . The tilded quantities are proportional to it; 4 r rf//2_4 , see Eq. (jCip . 

Note that although xi = % r ~ w 2 T d / 2 - 3 is one -loop order, the order parameter q(x) continues to be 

zero-loop order. 

We are now proceeding to compute the relevant nonanalytic terms to the free energy in (|14p : both 
and lnZg must be considered in this regime of dimensions. 

• Gaussian free energy: near infrared contribution. Eq. (|24|) remains valid, except replacing 112 by 
U2 in the second term, and exploiting the definition of U2 in (jCip : 



1 , 16 f d d p p 2 + l 27 4 _ o = 

In Z£ a = / — y — -f „ NO x r d / 2 w 4 «2 3 n 5 . 35 

niV G d ,/ (27r) d p 6 (p 2 + 2) 2 160d 2 v ; 



A detailed analysis based on Ref. show* t lint, unly highly subdoininant conl ribul ions 
mZ G 



-A? InZf? 1 ~ r d 1 — result from the far infrared regime p 2 ~ ~ ^h T<1 \ these are even smaller 



than the correction term from the near infrared regime which is ~ r d / 2 +i_ 

Nonanalytic temperature dependence occurs in of Eq. due to the tilded quartic couplings 
in the order parameter, see (|34p . For retrieving the relevant leading contributions, it is sufficient 
to take 

-^ (0) = ^ 2 r + ^ 3 r. 

thus neglecting terms which are smaller by U2/U2 ~ r 4 ~ d / 2 . A straightforward calculation provides 

_ J-tffD = C 2 . loop x w 2 r d - 5 + ±- w%\ 5 . (36) 
niv loL) 

(Note the lack of the term ~ T d l 2 .) The first term, which obviously dominates the subleading 
one ~ r d / 2+1 , consists of contributions like w^ l u 2 l q\. A consistent calculation of C2-I00P1 however, 
requires the 2-loop extension of the equation of state. Similarly, 2-loop corrections behind the 
Gaussian free energy in Eq. (|14p yield the same kind of term. 



16 



B. The replica symmetric case 

Similarly to the RSB case, we must include the one-loop term in Eq. (jlOp . and — keeping only the 
relevant terms — the RS equation of state reads: 

2rq + (n - 2)wq 2 + (n - 2)w-^ G 2 = 0, 

Q 



iV 
p 



with 



n(n - l)(n - 2)G 2 = 2(n - 2) 1 + (n - l)(n - 4) — -^L n(n - 3) ' 



P 2 + Al p 2 + Xa ' P 2 + Xr 

and for the RS masses see (|26p . Extracting the nonanalytic contribution from the third term is somewhat 
lengthy but straightforward; in the zero replica number limit it can be put into the form [see Eq. (IClj) ] 



p 



q 3 , n = 0. 



Note that the resultant RS equation of state is equivalent with the RSB one of Eq. (j33|) for x = Xq = x\. 

We are now in the position to analyse the two possible sources for the RS free energy having nonanalytic 
temperature dependence for 6 < d < 8: 

• As for the Gaussian part, Eq. (|27h remains valid for —jr \\xZq" even in this dimensional regime. 

• Using the RS equation of state from above, it is straightforward to see that the leading nonanalytic 
term to — ^y^ '* is proportional to r d ~ 3 . It has, however, a 2-loop character, and a consistent 
calculation would require extending to the next perturbative order. 



C. The free energy difference between the RSB ans RS phases 

Collecting pieces of information from previous sections and subsections, namely Eqs. (f22j) . (f27|) . ([35]) 
and (|36p. the leading terms which do not cancel in the free energy difference when 6 < d < 8 are the 
following: 

^RSB _ ^RS) = 

C 2 _ioop x w 2 r d ~ 3 + -'[ -^1 x r d / 2 n + -uT 2 rV + — (3 - d) w A u^n\ 

p dj (2n) d p 4 (p 2 + 2) 3 24 160d y ' 2 

[The notation C2_i 00 p was kept here for simplicity, although it may obviously differ from that defined in 

Eq. (]36p . For the definition of U2 see (ICip .] The line where the RSB and RS free energies coincide follows 

from this equation: 

n~w 2 r d / 2 - 3 , 6<d<8. 
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The proportionality constant is a one-loop integral, whose value can only be computed from the knowledge 
of the two-loop integral C2_ioop- 

VIII. THE STABILITY BOUNDARY OF THE RS PHASE 



It is well known from the famous finding of de Almeida and Thouless [10| that the mean field Ising 
spin glass in a homogeneous magnetic field enters the RSB phase along the boundary where the RS phase 
becomes unstable. This was later extended by Kondor [4 1 ] to the case when the magnetic field is zero, the 
replica number n, however, is finite: n essentially takes over the role of the magnetic field, and along an 
Almeida-Thouless (AT) line in the r — n plane, the RS phase becomes unstable. In [4| a simplified model 
- the "truncated" model with all the quartic couplings but u-i zero — was used, and it was shown that 
the RS and RSB free energies coincide along the instability line. 



In Ref. [12( the leading behaviour of the dangerous "replicon" mass Tr close to T c was expressed in 
terms of the exact cubic and quartic vertices yj exact and ti2 Xact , and also the exact order parameter g exact 

7 

as 

■p „ „ , ,exact exact „, exact exact 2 , c 

1 r = nw q — —11*2 q , d > o. 
o 

The zero replicon mode signals the instability of the RS phase providing for the AT-line: 

n AT = hw^^uf^q^. 

Substituting the bare values for the vertices and the mean field order parameter wq = r, agreement with 
the zero-loop result for the line of equal free energies of Eq. (|32p is found, thus reproducing Kondor's 
result. (Although the more generic model with all the quartic couplings is considered here.) 

Computing the one-loop correction of the replicon mass by standard perturbative methods is straight- 
forward for d > 8, although somewhat lengthy. Omitting the details, only the final result is displayed 
here: 



r R = {1 + [(n - l) Ul +u 3 + iu 4 ] I 4 + (n- 2)w 2 J 6 } n (wq e ™ ct ) + {- 



[n(n — 3)ui — 2u2 



+ - [(2n 3 - 5n 2 - 3n - 4)u\ - \2u\U2 + 4n(n - 3)uiu 3 + 24n(n - 3)«iu 4 - 18^ - 24n 2 n 3 - 48u 2 ^4 - 4u|] I 4 
+ i[-(2n 3 -6n 2 + 12n+16)ui+24u 2 -2(n 2 -8)n 3 +8n(n-6)n 4 ] w 2 I e -(2n 3 -9n 2 + 12n+16) w 4 I 8 \ q exact2 . 



7 The superscript "exact" is used to distinguish these quantities from the "bare" ones, which are the zeroth order contribu- 
tions to them. 
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This must be complemented by the equation of state, i.e. by the one-loop relationship between the order 
parameter and r: 

- 2r = (n - 2) j 1 + ~[(n - 2)m - 3u 2 - (n - 2)u 3 - (n 2 - n - 8)114] J4 - (n - 2)u> 2 / 6 1 (w<f xact ) 



+ 0(g e 



From these two equations, we can easily conclude the leading behaviour for n ~ t <C 1: 



r R 



1 + o + 3m 2 + ^3 + 4u 4 ) h - ^w 2 h 



TIT 



U2 H — (2u 2 + 10uiU2 + 15U 2 . 
3 



+ 8m 2 u 3 + 8u 2 u 4 + 2n|) J 4 + 8(iti - 2« 2 - u 3 )w 2 I G + 24w 4 I 8 



w- 2 r 2 . (37) 



The instability line is obtained from the condition = 0: 



L 9 

"-AT = ~[u2 + /d(A)] w~ r 



where the correction term has been defined in (|30j) . Comparing this expression with (132|) . we can conclude 
that the mean field type behaviour persists for d > 8: the RS phase becomes unstable where its free 
energy coincides with that of the RSB phase. (See, however, the discussion in the Conclusion part.) 



IX. CONCLUSION 



Two basic features of the mean field Ising spin glass were followed in this paper, by perturbatively 
taking into account the effect of the geometry of a high dimensional lattice (d certainly larger than 6). 
These properties are: 

• the anomalous sample to sample free energy fluctuations (considering only large deviations), 

• and the equality of the free energies of the replica symmetric and infinite step replica symmetry 
broken phases along the line (Almeida-Thouless line) where the RS phase becomes unstable. 

Both problems can be elaborated by studying the n-dependent free energy below the spin glass transition. 

As for the first item, perturbations break down the anomalous behaviour, and Gaussian large devi- 
ations take over the lead. As it was shown in Sec. IIVI Gaussian fluctuation is common for any ansatz 
of the order parameter with the property of replica eqivalence (not to be confused with replica symme- 
try), i.e. it must be a geometrical effect. As it was pointed out by G. Parisi [171 ]. Gaussian fluctuations 
always dominate whenever local interactions are inhomogeneous: this is certainly the case in the finite 
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dimensional geometry of a hypercubic lattice (but not for the SK model). The fact that locally inhomo- 



geneous interactions imply Gaussian large deviations of the free energy has been demonstrated recent 
on the Bethe lattice with finite connectivity, and continuously distributed quenched interactions 



y 



14- 



The anomalous n-dependence of the RSB free energy, however, does persist, although it is subleading in 
finite dimensions. 

The coincidence of the RS instability line and the line of equal free energies with the RSB phase is 
somewhat misterious even in the SK model: the free energy difference is of fifth order (in the double 
series in r and ra), as contrasted with the replicon mass wich is a second order quantity. We found that 
this feature of the mean field theory persists for d > 8, but in a rather nontrivial way (see the complicated 
correction term in Eq. (|30p ). We must emphasize that the two computations in Sees. IVT1 and IVnil are 
completely independent. This result gives an important support to the scenario, at least for d > 8, 
that replica symmetry must be broken by the infinite ultrametric hierarchy of the mean field spin glass 
proposed by Parisi. We must remember, however, that the condition of Eq. (|3ip was a priori assumed, 
and any conclusions depend on the validity of it. 

In the dimensional domain 6 < d < 8 the situation is more complicated. The AT line can be computed 
relatively easily, as it was in Ref. [131 ]: 

nAT = o w Ur z T ~ w T i 

O 

n 

see ()Cip and Eq. (47) of [13J]. Among the leading contributions to the free energy difference, which are 
now nonanalytic, there are terms which can be computed only by extending the calculation to two-loop 
order; as explained in Sec. IVIIi This seems to be unfeasible, and no a priori assumptions are available 
now. A real miracle would be the coincidence of the two lines, presumably due to several cancellations. 
One can speculate that otherwise the separation of the two lines might be explained by some non mean 
field scenarios: first order transition or replica symmetry breaking with replica equivalence, but not with 
the infinite step ultrametric structure. 
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APPENDIX A: TRACES OF POWERS OF THE GAUSSIAN MASS OPERATOR FOR A 
GENERIC RSB SCHEME WITH REPLICA EQUIVALENCE 



Traces appearing in the formula for \nZc, Eq. (116p . can be computed and arranged into the generic 
form 



1 



TrM fc ~ r k + 



expression of 



i2\aa ( i3\aa ( jA 



r E' 



.a/? 4 



These formulae are derived by two subsequent steps: 

• Substitution of the mass components from Q. 

• From the zero- loop equation of state H a p = 0, with m = — r in (jlip . and exploiting replica 
equivalence new, replica independent equations follow, like this: 



Ttf) m = -\w {4?) aa -\ 

2 o L 



In this way, we can get rid of terms where r and 4> are coupled. 

The following formulae are valid for any RSB scheme compatible with replica equivalence, not neces- 
sarily ultrametric: 

^Tr M = —(n— l)r \(2n - l) Ul + 3u 2 + (n + l)u 3 + in 2 - n + 4Wl (<f> 2 ) aa ; (Al) 

n 6 



^Tr M 2 = 2(n - l)r 2 + (n - 2)w 2 U 2 ) aa - -w\u x + 3u 2 + (n - 1)« 3 + (n 2 - n - 4Wl (<A 3 ) QQ 
n 3 

— u\ [{2n - 3)ui + 6u 2 + 2(n - l)it 3 + 2(n 2 - n - 8)u 4 ] (0 4 )" a H [2u 2 - 4(n - 2)mu 2 + 3u 2 

' ^ 18 
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2{n - 5)u 2 u 3 - 2(n 2 - n - 8)u 2 u 4 + 2uf\ ^ ^ + — [(3n - 4)u 2 + 12niu 2 + 8nm 3 - (n - 3) 



«3 



2(n 2 - n - 8)u 3 u 4 - n(n 2 - n - 8)u 2 4 ] ((0 2 ) Q °) 2 ; (A2) 



1 Tr M 3 = -4(n - l)r 3 + 2(n- 2)w 3 (c/> 3 ) QQ + ^w 2 [(2n - 7)m - 3u 2 - (n - l)u 3 - (n 2 - n + 4)1*4] (<t> 4 ) 



n 



+ w 2 [-2 Ul + {n + l)u 2 + 2n 3 ] ^ ^ + ^ [ ~ ( 6n ~ 19 M ~ 3n 2 + (n - 5)u 3 + 

(n 2 -n + 4H ((^ 2 n 2 + ...; 



1 



Tr M 4 = 8(n - l)r 4 + 



(3n - 4) + 8 £ + 2(n - 7) ((0 2 ) Qa ) 2 



+ .... 
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APPENDIX B: ZERO-LOOP ORDER PARAMETER OF THE r, u u u 2 , u 3 , u 4 MODEL 



For the model with higher than quartic invariants neglected — which is somewhat more general than 
the traditional "truncated" model where u\ = u% = 114 = — , the mean field equation of state follows 
from Eqs. (jlOp and (jlip by neglecting the 1-loop term and setting m = — r: 



ZD 



(Bl) 



(Replica equivalence was applied to get the last term.) The infinite-step ultrametric solution of this 
equation can be represented by the order parameter function q(x), see {i| and references therein, and x 
alls into the interval [n, 1] for n > small but finite. An application of the ultrametric replica algebra 
I9I ]. and extending the method used for the truncated model [20], a somewhat lengthy calculation leads 
to 



Z(*) 



K 

qi 



n < x < xq 
Xq < x < X\ 

%1 < x < 1, 



(B2) 



and the constant K is best expressed by x\ as 



1 + ^x\ 

U2 J- 



^ 1+^*1 



(B3) 



Continuity of q(x) at xq and x\ means that only two of the four parameters xq, x%, qo, and qi are 
independent. Choosing xq and q\ as such, the most economic way to give them as function of the model 
parameters is the inverse relationship: 

2 



n = o x o + 0(3), 



t = wqi 



1 

2 L 



1 



ui+u 2 - -(u 3 + nui) 



ql-\ U -^ql + 0^. 



(B4) 



w 



[Orders 0(. . . ) are understood in the sense of the double series in xq and q\. Notice the lack of 0(2) 
term in the first equation.] 

The invariants entering Eq. (|T71) can be computed by the methods presented in [19], and they have 
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the following expansions: 

(0 2 )«« = -q\ + ^ Xl q 2 + 0(5) = -w- 2 \t 2 + \w~ 2 {3 Ul + u 2 - u 3 - nu 4 )r 3 + . . 

{4> 3 ) aa = 2q\ - 2xiq\ + \x\q\ + — xfal + 0(6) = w~ 3 hr 3 + uT 2 (3ui + u 2 - u 3 - nu 4 )T 4 + ... 

+— w\ 2 3 n \ 



) 

E 

/3 



4\aa 



-3q$ + 4xigJ + 0(6) = -3-uT 4 r 4 + ..., 



uolB^ 4 , ^ 4 2 4 , , -4 4, 81 4 _ 4 5 

^ p = -9i + -aji?! - y^ <? + °( 6 ) = r + ^ w u 2 n > 



((^) fl,0 ) 2 = «f-|^ + 0(6) = u;- 4 r 4 + .... 



While deriving the above expansions, Eqs. (|B2p . ()B3P and (|B4p were frequently used. 



(B5) 



APPENDIX C: CALCULATION OF THE NONANALYTIC CONTRIBUTIONS FOR THE 

EQUATION OF STATE 

Extracting the leading nonanalytic contributions from the second term in Eq. (|1U|) becomes possible 
by the following simplifications: 



All the bare cubic vertices but W, 



w are negligible, see Eq. (jT2j) . 



We can use the bare propagators in the near infrared region 



15J ; this is equivalent by neglecting 



the bare quartic couplings u±, u 2 , 113, and U4 in the mass components in 



• The single replica sum arising, see below, is then approximated by evaluating it in the smallest 
ultrametric block of size x\ — 1, and x\ is obviously negligible here for small r and d > 6. 

We can write then 

x being the ultrametric distance between a and 0, and the parametrization of ultrametric matrices was 
used in the last formula 



lx 1 



Eq. (62) of Ref. 



15J] with y = x\ can now be used: 



2 (~iXX X 



{l + Qt + tt 2 -4i£)i a 



,(i +t) t\ 



(l + 2i) 2 * * (1 + 2£) 2 (l + 2t) 2 

where t = wqi/p 2 and t x = wq(x)/p 2 . Only the last two terms contribute to the nonanalytic temperature 
dependence we are looking for: the second term in Eq. (jlOp can then be written as 



1 



1 



ui ~ ) Qiqyx) + -u 2 q{x) 
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with 



Ul 



oo 

1= 1R f d d p p 2 + l 4 d/2 _ 4 

no 

d d p 1 



(CI) 



u 2 = 24 



J (2vr) d p 4 (p 2 + 2) : 



4 _d/2-4 



■ to r 
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